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CMB Distortions from Damping of Acoustic Waves Produced by Cosmic Strings 

Hiroyuki TashiroQ Eray Sabancilarlll and Tanmay Vachaspat^l 

Physics Department, Arizona State University, Tempe, Arizona 85287, USA. 

We study diffusion damping of acoustic waves in the photon-baryon fluid due to cosmic strings, 
and calculate the induced n- and y-type spectral distortions of the cosmic microwave background. 
For cosmic strings with tension within current bounds, their contribution to the spectral distortions 
is subdominant compared to the distortions from primordial density perturbations. 
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Spectral distortions of the cosmic microwave back- 
ground (CMB) are an important tool to study the ther- 
mal history of the universe at epochs earlier than re- 
combination [H, m. COBE-FIRAS measurements of the 
spectrum of the CMB show a blackbody spectrum to high 
accuracy, and result in upper bounds on the spectral dis- 
tortion parameters, < 9 x 10^^ and y< 1.5 x 10"^ 
[3, 01 • Proposed missions such as PIXIE ^ , will be able 
to detect distortion parameters, and y, up to levels of 
10~^ with 5-f7 significance, spurring theoretical investi- 
gation of early universe models and their predictions for 
spectral distortions 

There are various physical mechanisms that can in- 
duce spectral distortions, including dissipation of acous- 
tic waves in the photon-baryon fluid created by the pri- 
mordial density perturbations via Silk damping 14|-|l6|. 
decaying particles 17|, electromagnetic radiation from 
superconducting strings Q , Hawking radiation from pri- 



mordial black holes [18[ and diffusion of magnetic fields 
[l^ . Damping of acoustic waves produced by primor- 
dial density perturbations have been studied extensively 
0i @]- In this paper, we study cosmic strings as an 
additional source of acoustic waves, and calculate the re- 
sulting CMB distortions due to the diffusion damping of 
these acoustic waves. 

The reason cosmic strings induce CMB spectral dis- 
tortions is that they are relativistic, massive objects, and 
they stir up the primordial fluid by their time dependent 
gravitational potential to generate acoustic waves. This 
is a continuous process from the time that the cosmic 
strings were produced. Acoustic waves dissipated due to 
Silk damping in the redshift range 10^ > z > 5x 10''^ 
produce /i— distortions of the CMB; waves dissipated in 
the range 5 x 10^ > z > 1100 produce y— distortions. 

The production and evolution of acoustic waves is ob- 
tained by solving the continuity and Euler equations 
for the longitudinal (scalar) modes for the photon and 
baryon fluids in Fourier space. We follow the treatment 
in Ref. [20| with the inclusion of a source term that de- 
scribes the effect of strings. The equations for comoving 



mode number k are 
— ——kv^, 
k, k 

Sb = -kvb, 



t{v^ ~ Vb) + S{k,ri), 
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Vb) + S{k,r)). 



(1) 

(2) 
(3) 
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The subscripts 7 and b correspond to photon and baryon 
fluids. Pi is the background energy density {i — 7,6), 5i 
is the density contrast, Vi = k ■ Vi and ni = k ■ iZi are the 
longitudinal modes of the velocity and the anisotropic 
stress, S{k,r]) = k ■ S{k,ri) is the longitudinal mode of 
the source term, a is the scale factor, the overdot rep- 
resents differentiation with respect to conformal time 77, 
f = aneCT with the Thomson scattering cross section, 
(Tt, and the electron number density, rig. We are inter- 
ested in the cosmic string contribution to spectral distor- 
tions, and so we will use vanishing density perturbations 
and fluid velocities as initial conditions. 

The photon and baryon fluids are tightly coupled on 
the scales of interest since k/f <^ 1. Then we expand the 
difference in velocities to second order in fc/r, 

v^-Vb^ (fc/f)/(T) + (k/rfgir), (5) 

where /(t) and g{T) are functions that we will determine. 
We will also ignore terms involving a /a since the damping 
process is much faster than the Hubble expansion. 

To first order, Vb ^ - {k/f)f, and Eqs. ^ and (g]) 
give 

^7 = 7^7 - 



4 
k 



S{k,v), 



= -^f + ^(fc,??), 

where R = 3pb/'ip^ oc a. To lowest order in 
drop the anisotropic stress term because 

96 fc 



(6) 
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Vb, we 
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in the tight coupling approximation. Eliminating / from 
both equations, we obtain 

k 

— 



4(1 



Sik,r^). 



(9) 
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Comparing the result with Eq. ([6]), we get 

R 



4(1 + 



According to Eqs. ([T]) and ([9]), we have 



3(1 + 



^7 - -kS{k,T]), 



(10) 



(11) 



which shows that the source term will generate acoustic 
oscillations that propagate at the sound velocity = 
1/^3(1 + 

To next order, we include the anisotropic stress term 
([8]) and obtain 



fc „ 16 P 

4 



^ 4 ^ 45 T ^ 
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(^7 + -ff 



+ S, (12) 
+ 5. (13) 



The dilference of these equations yields g, and then, using 
Eq. (dJ, we get 



k 3 
i'7 = 4(i + i^) '^'y + 2^'^'^'^ ^ 



where 



1 



16 



1 + i? 15 



67-(l + R) 

Using Eqs. ([Ij and (ITi|) . we finally obtain 
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^u,^2eD5^~-kS{k,r^). 



(14) 



(15) 



(16) 



An approximate analytical solution of Eq. (jl6p is 



-k^'/kl 



cos(fcr,) / S'(fe,ry')sin(/cr;) 



sin(fcrs) / S{k,T]') cos{kr'g) 



, (17) 



where the sound horizon r^, and the Silk damping scale 
l/ko t21i] are 



r's = rs{rii,r)') = drj Cs, 

Jr)i 

■n' 

D ^ 



/c-2 ^ / dr^ D 



(18) 
(19) 



Vi 



Note that the solution given in Eq. pT|) satisfies Eq. (fT6|) 
in the tight coupling approximation i.e., we have ignored 
terms of order (fc/f)^ and higher. We have also ignored 
terms containing Cg, R and H = a/ a, since the rate of 
damping of acoustic waves is a lot faster than the Hubble 
expansion rate. The expression in Eq. ()17p applies to a 
single source and, by linear superposition, to a collection 
of individual sources. 



The energy density of acoustic ( "sound" ) waves is given 

by BEl 

= ^p(l + 3c^)(|<5,(a.,^)p) « ^p{\S,ix,rj)\'), (20) 

where we approximate w 1/3. The averaged density 
fluctuation over a volume V ^ oo on the right-hand side 
of Eq. (|20p is defined in terms of Fourier modes by 



cPk 



V J (27r)3 



(21) 



The rate at which energy density is dissipated from 
acoustic waves due to Silk damping is found by substi- 
tuting Eq. (fT7|) in Eq. (|20|) . and then by taking the time 
derivative and only retaining the contribution due to Silk 
damping, i.e., only keeping the term coming from the ex- 
ponential damping factor as 



d^k 
(2^ 



k" \5^{k.v)? 



(22) 



where 5{k, rj) is given by Eq. (jl7p and S{k, rf) is the con- 
tribution from all sources. 

The spectral p- and y-distortions from these sources 
are given by 11 1 



l)fr< 



P = 1A I dr] , 



1 



drj 



Qiv) 

3p ■ 



(23) 
(24) 



where ij^, ?7frcczc and rj^^c correspond to epochs 
2 X lO'', 2;fj-ccz 



5 X 10 and z^, 



1100. Note that 



the factor 1/3 in Eqs. (1^51) and (1^^ comes from the fact 
that only 1/3 of the dissipated energy creates the CMB 
distortions while the remaining 2/3 is used to raise the 
average CMB temperature, which cannot be observed di- 
rectly 0- 

We now consider distortions due to long strings. We 
will adopt the scaling solution for the cosmic string net- 
work suggested by Nambu-Goto simulations The 
number density of long strings (in comoving volume) is 



nio 



K a 



(25) 



where a is the scale factor, n = 11 for the radiation domi- 
nated era, and k = 15 for the matter dominated era from 



the recent simulation results (22| . The dependence of the 
density on reconnection probability p is expected to be of 
the above form [23, 24[. The reconnection probability is 
unity for ordinary gauge strings, and can be smaller than 
unity for cosmic F- and D-strings, e.g., 10~^ ^ P ^ 1 j25| . 

Numerical simulations also show that long strings are 
not straight, but have wiggles on them. The difference of 
a wiggly string from a straight one is that its tension, T, 



3 



is not equal to its mass per unit length, jls ■ They satisfy 
the equation of state given by Tjis = A*si where Hs is the 
bare mass per unit length of the string. The gravitational 
force due to a wiggly long string is F cx {p.^ — T)/r at a 
distance r from the string. Thus, for a single long wiggly 
string lying along the t/-axis and moving along the z- 
axis with comoving velocity Vg, the gravitational force 
due to the string at the point x = (x,y,z) in comoving 



coordinates is [26|, |27 1 



{x - Vs-qz) 



\x 



(26) 



where, w = (fig ~ T)/ ^.g is the parameter that controls 
the wiggleness, fj = — rji is the elapsed comoving time, 
z is the unit vector along the 2:-axis, and the superscript 
denotes that this is the source term for a single long string 
with label i. The location, orientation and velocity of the 
string will be different for different long strings. We will 
average over the locations and orientations, and use the 
root-mean-square velocity as seen by simulations. 

The Fourier component of the longitudinal mode of the 
source term, Eq. (pS)) . is 



'longV 



{k,r,)^ I d'xk-S^,i^{x,Tj)e 



ik-x 



= S2i7:^Gfis e™^'?'^^- [1 - Jo(fc,77)] , (27) 

where k'^ = + k'^., ~ + k^, the subscript D in So 
refers to the Dirac delta function, and we assume that the 
gravitational force due to the long string does not prop- 
agate beyond the horizon size rj. The total source due to 
all long strings will be given by linear superposition 



5iong(fc,ry)-5]5^,^;;g(fc,,7). 



(28) 



The contributions from different long strings will add in- 
coherently, and Eq. (1^ gives 



(29) 

where niong is given by Eq. ([25]) and J^^"* is the density 
fluctuation [see Eq. ()17p for a single source 

We now numerically evaluate the integrals in Eqs. 
(|24| with Q from ^ and S'^^^ from ^ to get 



1. X 10 



-13 



p \10-' 



2. X 10"'* — 



p \10 



(30) 



(31) 



The COBE-FIRAS measurements of the CMB spectrum 
constrain |^| < 9 x 10^ and |?;| < 1.5 x 10"^ [3], and the 
detection limits by PIXIE are ^ > 5 x 10~^ and y > 10~^ 
Therefore, the distortions due to long cosmic strings 



are too small to be detected at PIXIE for strings with 
tensions below current bounds, G/ig < 10"^. 

Having discussed CMB spectral distortions from long 
strings, we now turn to cosmic string loops. For a loop 
of size L and tension G/jg moving along the z-axis with 
comoving velocity Vg, the gravitational force at the point 
X = (a;, y, z) in comoving coordinates is 



*^loop 



G^gL [X - VsT]z) 



\x - Vgr]z\ 



(32) 



Here a is the scale factor at conformal time 77, fj = r] — rji 
is the elapsed comoving time, rji is the conformal time at 
loop production, L is the physical loop size. Note that 
we model the loop as a point source. This is a good ap- 
proximation at distances larger than the loop size, where 
most of the acoustic wave energy resides. 

The Fourier transform of the longitudinal component 
of the source term (l32l) is 



(33) 



S'ioop(fc, V,L) = J d^x k ■ 5'ioop(a;, 77, L) e''' '- 



sin(fc?7) 



k^rj 



where k^ = kl + ky + k'^. 

The formula for loops analogous to for long strings 
can be written as 



Qioop(?7) = Y 



d^k ., 



{2n) 



(34) 



where dnioops is the number density of loops of length 
between L and L + dL. Note that S^{k,ri) cx L and so 
the last factor in ([M)) is independent of L. 

The number density of loops has not been completely 
resolve43- However, there is substantial evidence for a 
scaling loop distribution in recent simulations (22j . In 
what follows, we adopt the results of the Nambu-Goto 
simulation results from Ref. , where it has been found 
that a significant fraction of the energy density that 
goes into loop production peaks at loops of initial size 
Li ~ ati, where a ~ 0.1 and U is the cosmic time at 
formation. There is a distribution of loop velocities with 
an average of Vi ~ 0.3 from the same simulations. The 
number density of cosmic string loops of length i as a 
function of time t (in comoving volume) is 28] 



dL 



dnioopiL,t) = 



I ty^{L + TGfigtf/^ ' 



t < t 



cqj 



CLdL 



f^iL + TGiigt) 



2 ' 



t > ieq, (35) 



^ See Ref. [22| and the references therein for a detailed discussion 
of the subject. 
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2 [22], p ^ 1 is the reconnection probability, 

and 1 for L > 



TGfist) for L < ate 



where kl 
Cl is \/tcq/{L 

Note that the distribution given by Eq. (|35l) takes into 
account the energy fosses of foops due to gravitational 
radiation, i.e., as a loop emits gravitational waves, its 
length decreases as L{t) « L{ti) — TG^gt, where its life- 
time is - L/{TGiJLs) with T ~ 50. 

The loop distribution in Eq. ([55]) yields 



dnioopL'^ 



P 



.3/2 



tcq/t^ + Va/t , 



i>tcq, (36) 



and numerical evaluation of the last term in (1341) finally 
gives us 



—--"KM 



(37) 
(38) 



Therefore, the string loops cannot produce detectable 
CMB distortion either. 

In this letter, we have investigated the evolution of 
acoustic waves produced by long cosmic strings and loops 
in the photon-baryon fluid. Both long strings and loops 
can create acoustic waves gravitationally, and the ampli- 
tude is proportional to the string tension, G/is • We have 
also evaluated the CMB distortions due to the dissipa- 
tion of the acoustic waves. However, primordial pertur- 
bations also produce acoustic waves and their dissipation 
is guaranteed to produce CMB distortion. The acoustic 
waves from primordial perturbations, with WMAP cos- 
mological parameters, can produce distortions of order, 
II ~ 10~* and y ~ 10~^. Therefore the distortions due to 
long cosmic strings are completely overwhelmed by those 
produced by the primordial acoustic waves. On the other 
hand, long cosmic strings with Gfj, ~ 10~^ can produce 
CMB anisotropics with amplitude roughly one hundredth 
of that due to primordial acoustic waves at the scale 
of the sound horizon at the recombination epoch p9| . 
Considering this fact with Eq. ([3T|) . the obtained CMB 
distortions due to long strings seem to be much smaller 
than what we might have naively expected. The rea- 
son is that the amplitude of the acoustic waves produced 
by long strings have a highly red spectrum (spectral in- 
dex is equal to —4) below the cosmological horizon scale. 
Therefore, the amplitude of the acoustic waves around 
the Silk scale is much smaller than that of the primordial 
acoustic waves. We find that the acoustic wave ampli- 
tude due to cosmic string loops is peaked at the scale 
1/A; ~ at/ a. The amplitude is roughly proportional to 
k"^ below the peak scale. Thus, the acoustic waves do not 
have large enough amplitude around the Silk scale to cre- 
ate detectable CMB distortions. However, we would like 
to note that if cosmic strings are superconducting, they 



can produce detectable CMB distortions [8|. Therefore, 
CMB distortions can still be useful for studying the ef- 
fects of cosmic strings. 

We would like to emphasize that we only considered 
fixed loop velocity, Vs ~ 0.3, which is below the sound 
spee d, Cs . In reality there is a distribution of loop veloc- 
ities |30l . l3l| with smaller loops being born with larger 
velocities |32| . When Vs > Cg, our analysis does not 
hold, since Eq. ^TE\\ shows that there are values of the 
wavenumber k at which resonance can occur. In addi- 
tion, supersonic strings and loops can produce shocks 
|27. . ,33 ] and these are not included in our analysis. How- 
ever, since cosmic expansion decreases the loop velocity, 
small loops are expected to be supersonic only for a finite 
period of time. Further, small fast-moving loops carry 
only a small fraction of the total energy in loops. So we 
do not expect supersonic effects to change our estimates 
significantly. 

We would like to thank Eiichiro Komatsu and Matias 
Zaldarriaga for useful discussions. This work was sup- 
ported in part by the Department of Energy, National 
Science Foundation Grant No. PHY-0854827 and Cos- 
mology Initiative at Arizona State University. 
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